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ABSTRACT 


Most  studies  of  upper  ocean  circulation,  tropospheric  dynamics,  and  air-sea  interaction  require 
the  measurement  of  near-surface  vector  winds  over  the  oceans  under  all-weather  conditions. 
Radar  scatterometers  may  be  used  to  find  the  surface  wind  vector.  The  scattering  coefficient  of 
the  ocean  usually  is  assumed  to  take  the  form 

o°  =  AuT 


where  u  is  the  wind  speed. 


To  determine  the  correct  wind  exponent  y  by  microwave  scatterometers  is  therefore  very 
important.  Many  experiments  on  scattering  from  the  oceans  provide  similar  but  different  values 
of  y. 

Most  experimenters  apply  the  ordinary  least-square  regression  (OLR)  to  extract  y  from 
measurements.  However,  this  technique  assumes  that  only  one  variable  may  be  in  error,  so  such 
a  result  may  not  be  statistically  correct.  As  an  alternative,  orthogonal  regression  (OR)  provides 
the  best  results  when  both  variables  have  similar  scales,  the  wind  exponent  y  is  discussed  in  two 
cases  with  both  OR  and  OLR.  Case  one  is  when  both  wind  direction  and  wave  direction  are 
same,  in  the  Up  direction,  and  case  two  when  the  two  direction  are  orthogonal.  For  case  one, 
OR  and  OLR  were  found  to  have  close  results;  however  results  are  quite  different  for  the  second 


It  is  necessary  to  find  optimum  values  of  y  at  different  incident  angles  and  dirctions.  The 
optimum  y  can  be  used  to  obtain  the  ocean  surface  winds  peed.  The  bilinear  regression  described 
in  the  report  is  a  simple  but  effective  way  to  show  he  relation  between  y  and  incidence  angles. 
Results  using  OR  provided  a  very  good  fit  to  the  bilinear  regression  analysis  of  various 
experimental  data. 
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CHAPTER  ONE 


INTRODUCTION 


1-1.  History  of  experiments  for  radar  ocean  wind  measurements  (1968--). 

The  concept  of  inverting  radar  scattering  amplitudes  to  determine  the  speed  of  the  ocean  surface 
winds  using  a  scatterometer  was  proposed  by  Moore  and  Pierson  in  1966  [l](Moore,  R.K, 
WJ.PiersonJr.  1966).  Since  then  many  experiments  have  been  conducted  to  determine  the  radar 
response  under  varying  wave  and  wind  conditions,  and  the  theory  of  radar  return  from  the  sea 
has  advanced  significantly. 

One  of  the  most  extensive  experimental  programs  was  initiated  during  the  1950s  by  the  U.S. 
Naval  Research  Laboratory.  This  program  conducted  experiments  using  a  multi-frequency 
airborne  radar  during  a  variety  of  sea  conditions  and  at  various  radar  look  angles.  The  results 
of  the  1950s  experiments  were  reported  in  [2] (MacDonald,  1956). 

During  the  1960s,  the  Naval  Research  Laboratory  (NRL)  continued  the  extensive  aircraft 
measurements  program  to  provide  quantitative  information  on  the  parametric  behavior  of  the 
electromagnetic  scattering  coefficient  (RCS)  o°  of  the  ocean.  Measurements  were  obtained  as  a 
function  of  polarization,  incidence  angle,  and  azimuth  angle  using  4-frequency  pulse  radars 
operating  at  0.4, 1.2,  3(1950s),  4.4(1960s)  and  9.0  GHz.  Other  aircraft  measurements  were  later 
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performed  by  NASA,  New  York  University  and  The  University  of  Kansas  using  a  NASA 
Johnson  Space  Center  radar  scatterometer  at  13.3  GHz  (  a  fan-beam  Doppler  radar).  Analysis 
of  NASA  observations  at  13.3  GHz  and  NRL  observations  at  8.9  GHz  and  4.4  GHz  show  the 
variability  with  wind  speed.  A  detailed  analysis  is  given  in  [3}(Claassen,  Moore.  1972).  The 
NASA-JSC  data  sets  were  used  to  propose  a  microwave  scatterometer  technique  for  remotely 
sensing  surface  windspeed  over  the  ocean. 

More  detailed  correspondence  between  microwave  scatterometer  backacatter  crosa-aection,  o°,  and 
measured  surface  wind  was  established  in  the  1970s  using  an  airborne  scatterometer  flying  near 
buoys  and  ships.  This  was  extended  to  satellite  data  when  NASA  launched  Skylab  in  1973  and 
SEASAT  in  1978.  During  this  decade,  the  main  experiments  were: 

*  AAFE  program:  An  aircraft  microwave  scatterometer  was  developed  at  NASA  Langley 
Research  Center  under  the  Advanced  Applications  Flight  Experiments  (AAFE)  program.  These 
experiments  evaluated  the  viability  of  radar  remote  sensing  of  the  surface  wind  vector(both  speed 
and  direction)  using  a  circle  flight  experimental  technique  suggested  by  Jones.  The  results  are 
shown  in  Fig.1-1,  and  1-2  [4](Ulaby,  F.T.,  Moore,R.K.,Fung,A.K.,1986). 

♦Skylab  Measurements:  SKYLAB.  launched  in  1973,  provided  the  first  opportunity  to  observe 
the  oceans  from  space.  The  SKYLAB  S-193  instrument  contained  a  radar  scatterometer  / 
radiometer  at  13.9  GHz  and  a  radar  altimeter  at  the  same  frequency.  With  the  great  speed  and 
vast  coverage  of  a  space  borne  instrument,  only  a  few  days  observations  were  more  than  all  the 
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(a) 


Figure  1-1  Example  of  circle-flight  observation*  of  the  azimuthal  variation  of  o*  at  13. 
GHz,  along  with  a  regression  line  of  the  form  of  Eq.  11,46  (from  Moore,  el  al.,  1978).  <t>  •  < 
corresponds  to  the  upwind  direction. 


Scotterlng  Coefficient  a*  <dB> 


Figure  1-2  Comparison  of  AAFE  Radscat  J3.9  GHz  <r*  measurements  with  theory  of 
Chan  and  Fung  (1977)  (frdm  Moore,  el  o!.,  1978). 


preceding  programs  combined  [5]-[10].  However,  since  the  spacecraft  moves  so  fast  and  the 
time  available  for  this  particular  experiment  in  the  complex  of  SKYLAB  experiments  was 
limited,  very  few  points  were  obtained  where  "surface-truth"  wind  measurements  were  available 
from  the  surface.  Comparison  of  Skylab  data  with  an  underflight  by  the  AAFE  Radscat 
instrument  showed  they  were  within  about  0.5  dB  of  each  other.  Comparison  was  also  made 
with  objective  analyses  of  reports  sent  in  by  900  ships  on  the  North  Atlantic.  The  ship 
measurements  are  not  very  reliable,  but  the  results  at  least  show  that  the  value  of  0°  increases 
with  wind  speed  at  incidence  angles  larger  than  10  degrees,  but  decreases  at  vertical  incidence. 
See  Fig.1-3  and  1-4.  [4]  (R.K.Moore  vol.iii)) 

♦Seasat  Measurements  of  Ocean  Backscattering:  The  SEASAT-A  satellite  was  launched  in  June 
1978.  The  goal  was  to  develop  meteorological  and  oceanographic  measurements  from  space  with 
microwaves.  The  measurement  lasted  only  99  days,  but  it  provided  reliable  scatterometer  data 
and  SAR  imagery  under  a  variety  of  wind  and  wave  conditions.  To  compare  measurements  made 
by  SEASAT  sensors  with  contemporaneous  surface  observations  from  a  number  of  research 
vessels  and  moored  buoys,  two  major  experiments  involving  extensive  surface  measurements 
were  conducted:  the  GOASEX  experiment  in  the  Gulf  of  Alaska,  and  the  JASIN  experiment  in 
the  North  Atlantic  between  the  British  Isles  and  Iceland.  The  details  of  the  two  campaigns  were 
presented  in  special  issues  of  JGR,  1982;  and  JGR,  1983.  The  results  show  that  the 
measurements  of  wind  exponents  from  SEASAT  are  comparable  with  the  previous  aircraft 
experiments  and  Skylab  campaign  [Fig.1-5,1-6,  and  1-71  [4](RKM,  Vol  iii)]. 


Figure  1-3 (a)  Regression  fit  for  Skylab  tr*  measurements  at  32*  angle  of  incidence  (from 
Moore  and  Young,  1977). 


20 


0  5  10  15 

Scattering  Coefficient  0°  (dB) 


rigure  1-3  (b)  Regression  fit  for  Skylab  o®  measurements  at  1*  angle  of  incidence. 


(from  Moore  and  Young,  1977). 


Kind  Speed  (ms-1) 

Figure  V-4  Tower-based  measurements  of  the  wind  speed  response  of  scattering  at  10 
GHz.  Notelhat  y  -  ct'VcosS  is  recorded,  and  that  the  angles  given  are  the  grazing  angles, 
j.e.,  the  complements  of  the  incidence  angles  (from  dc  Loor  and  Hoogeboom,  1983). 


In  the  1980s  no  space  borne  scatteromenter  flew,  but  there  were  many  airborne  campaigns.  In 
this  period  the  major  experiments  were: 

TOSCANE  campaign  by  IFREMER  began  in  1983:  The  purpose  was  to  study  the  structure  of 
ocean  surface  wind  and  to  develop  an  array  of  instrumented  buoys  to  be  deployed  at  sea  during 
the  post  calibration  phase  of  the  ERS-1  C-band  wind  scatterometer.  The  Erst  campaign  named 
PROMESS/TOSCANE I  was  conducted  over  the  Atlantic  ocean  close  to  Brest,  France,  in  1984. 
It  dealt  with  ship  wind  measurements  and  with  airborne  scatterometer  calibration.  In  1985  (from 
Feb.  14  to  April  17),  a  second  campaign  called  TOSCANE  was  performed  near  southern  Brittany 
in  France.  The  goal  of  the  campaign  was  to  compare  the  aircraft  scatterometer  data  to  data  from 
an  array  of  well-calibrated  buoys  for  the  ERS-1  satellite  wind  scatterometer  calibration.  To 
estimate  the  performance  of  the  planned  ERS-1  C-  band  W  polarization  wind  scatterometer, 
several  airborne  scatterometer  measurements  were  carried  out  in  the  ESA’S  Earth  observation 
program  with  the  multifrequency  air-bome  scatterometer  DUTSCAT  in  the  TOSCANE  T 
campaign  in  1987.  Long  developed  a  scattering  model,  CMOD1,  using  the  DUTSCAT  signature 
of  the  ocean  database,  [11]  (Unal,  C.M.H.  et  al.,1991)  and  [12]  (Daniault,  N.,1988). 

FASINEX:  From  winter  of  1985  to  spring  of  1986,  the  Frontal  Air-Sea  Interaction  Experiment 
(FASINEX)  was  conducted  in  the  subtropical  convergence  zone  southwest  of  Bermuda.  The  goal 
of  the  experiment  was  more  understanding  of  the  relationships  between  radar  backscatter  from 
the  ocean  surface  and  the  parameters  of  geophysical  theory,  which  include  the  ocean  wind  stress, 
air-sea  temperature  difference,  and  sea-surface  wave  spectrum.  In  the  experiment,  the 
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observations  and  measurements  were  performed  by  buoys,  ships,  aircraft  and  spacecraft.  Most 
of  the  data  of  the  NASA  Jet  Propulsion  Laboratory  Active  Microwave  Scatterometer  (AMSCAT) 
came  from  10  measurement  flights  by  the  Cl  30  aircraft.  The  AMSCAT  operated  at  Ku  band. 
For  each  measurement,  the  aircraft  traveled  12km  along  a  straight  flight  path  at  1100  m  altitude. 
The  radar  incidence  angle  was  from  0  to  60  deg.  and  the  azimuthal  scan  was  about  330  deg. 
The  ocean  surface  wind  varied  from  2  to  20  m/s.  The  total  time  of  data  collected  was  30  hours 
[13]  (Weissman.  D.E,  1990)  and  [14](Plant,  WJ.,1990) 

NORCSEX:  In  March  1988.  an  ERS-1  prelaunch  experiment  on  the  Norwegian  Continental  Shelf 
(NORCSEX)  was  conducted.  The  basic  objective  of  the  experiment  was  to  investigate  the 
imaging  capability  of  SAR  on  the  ocean  current  surface.  A  secondary  objective  of  NORCSEX 
was  to  assess  the  capability  of  a  C-band  SAR  to  measure  ocean  surface  wind  and  waves.  As  the 
most  critical  parameter  for  remote  sensing  (SAR  and  scatterometer).  the  wind  stress  was 
measured  by  buoys  and  ships  during  NORCSEX.[15],[16],[17]. 

TOWARD:  The  Tower  Ocean  Wave  and  Radar  Dependence  Experiment  was  conducted  from 
October  1984  to  January  1986  at  the  Naval  Ocean  Systems  Center  Tower,  located  offshore  of 
Mission  Beach,  San  Diego,  California.  The  general  goal  of  the  experiment  was  to  provide  an 
experimental  basis  to  resolve  the  conflicting  theories  on  synthetic-aperture  radar  (SAR)  imaging 
of  the  ocean  surface.  The  specific  objectives  stated  for  the  TOWARD  experiment  included  (1) 
investigation  of  the  hydrodynamics  of  short  waves  and  their  modulation  by  long  waves,  (2) 
assessment  of  the  assumptions  stipulated  in  radar  backscatter  theory  that  are  used  in  SAR 
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Figure  1-6  Downwind  speed  response  exponents  from 
Offerer.;  experiments. 
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different  experiments. 


imaging  of  the  ocean  surface,  and  (3)  developing  a  verifiable  theory  for  SAJI  imaging  of  the 
ocean  surface  [18](Shemdin,O.H.,1988). 

SAXON -CLT:  A  significant  achievement  from  TOWARD  is  the  determination  that  none  of  the 
then  available  theories  on  SAR  imaging  of  long  surface  waves  could  explain  all  the  SAR 
observations  satisfactorily,  particularly  at  higher  frequencies  like  X  band.  To  improve  the 
theoretical  models,  the  SAXON-CLT  experiment  was  executed  during  September  to  October  1988 
at  the  Chesapeake  Light  Tower  (CLT),  located  15  miles  east  of  Virginia  Beach,  VA.  The 
primary  objectives  were:  (1)  investigate  the  hydrodynamic  modulation  of  short  gravity  and 
capillary  waves  by  long  surface  waves  and  internal  waves,  (2)  determine  experimentally  the 
modulations  in  radar  waves,  (3)  test  SAR  imaging  theories  at  X-Band  and  higher  frequencies,  and 
(4)  investigate  the  influence  of  the  microlayer  on  both  active  and  passive  remote  sensing 
techniques  [19](Popstefanja  and  McIntosh,  1989). 

Tiie  SAXON-FPN  experiment  was  conducted  in  November  of  1990  with  similar  objectives. 

1-2.  Background  of  SAXON  experiment 

The  Synthetic  Aperture  X-Band  Ocean  Nonlinear  Experiment  (SAXON)  was  by  American  and 
German  experimenters  in  the  North  Sea  on  and  near  the  German  Forschungsplatform  Nordsee 
(FPN)  in  November  of  1990.  The  primary  objective  of  the  experiment  was  to  investigate  the 
influence  of  the  nonlinear  parameters  of  the  sea  surface  on  the  radar  backscatter  and  the  SAR 
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Figure  1-8  (a)  The  range  of  windspeed  in  SAXON- FPN 


imagery  at  a  variety  of  microwave  frequencies  under  high  sea-state  conditions. 

The  Radar  Systems  and  Remote  Sensing  Laboratory  (RSL)  of  the  University  of  Kansas  joined 
this  experiment  by  using  a  switched-beam  vector-slope  gauge/scatterometer(VSG)  (Ka-band)  and 
a  dual-frequency  X-  and  C-band  FM-CW  radar.  One  objective  was  to  improve  knowledge  of  the 
ocean  -surface  wind  exponent  7,  which  is  given  by  the  empirical  model, 

o°  =  A  ut  (2) 

where  a0  is  the  radar  scattering  cross-section,  u  is  the  ocean  surface  windspeed  and  A  is  a 
constant.  Figure  1-8  (a)  and  (b)  show  the  range  of  wind  conditions  encountered  during  the 
month. 

In  Chapter  2  of  this  report,  the  principles  of  windspeed  exponent  measurement  are  reviewed. 
Because  of  significant  errors  in  windspeed  measurement,  and  possibly  in  radar  measurement,  care 
in  using  regression  is  necessary.  Two  algorithms,  orthogonal  regression  and  ordinary  linear 
regression,  are  studied.  The  measurements  errors  are  also  presented. 

In  Chapter  3,  as  an  example  of  the  windspeed  exponent  at  35  GHz.  the  SAXON  experiments  are 
studied  in  more  detail.  In  this  chapter,  power  measurement  and  data  classification  are  discussed 
and  the  results  from  orthogonal  regression  and  ordinary  linear  regression  are  compared. 
Chapter  4  discusses  variation  of  wind  exponent  with  angle  of  incidence  in  detail.  Bilinear 
regression  is  derived  and  used  to  explain  the  data.  Finally,  Chapter  5  states  the  conclusions. 
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CHAPTER  TWO 


W1NDSPEED  EXPONENT  MEASUREMENT 

The  data  from  many  aircraft  scatterometer  and  spacebome  scatteometer  experiments  demonstrated 
that  the  radar  backscattering  coefficient  o°  over  the  ocean  can  be  used  to  infer  near-surface 
oceanic  winds  through  a  geophysical  model  function.  In  our  case,  we  use  the  commonly 
accepted  empirical  model  given  by 

o°  »  AuY  (3) 

Equation  (2)  can  be  expressed  in  log  arthmic  form  by 

logo0  «  log  A  +  ylogu  H) 

One  can  see  that  for  this  model,  log  o°  is  a  linear  function  of  log  u.  Therefore,  to  determine  the 
coefficient  y,  the  slope  of  the  linear  equation  (3),  is  the  most  important. 

2-1.  Linear  Regression  Methods 

We  are  concerned  with  the  analysis  of  the  relationship  between  two  variables,  x  and  y,  where  one 
variable  (y)  is  statistically  dependent  on  the  other  (x).  Both  variables  are  measured  on  a  fixed 
scale,  and  the  relationship  between  them  is  assumed  to  be  linear.  A  me.  iod  to  fit  a  first-degree 
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polynomial  function  ( the  linear  function  Y  =  b  X  +  b,  is  considered  here  )  to  a  set  of  observed 
paired  scores  is  required.  In  our  case,  the  Y  variable  is  <f  in  dB  and  the  X  variable  is  log  of 
ocean  wind  speed  u  ,  incidence  angle  or  some  other  geophysical  quantity. 

Most  experimenters  use  ordinary  linear  regression  (OLR)  of  the  radar  scattering  cross-section  a° 
vs  ocean  surface  windspeed  u.  The  critical  assumption  of  the  OLR  method  is  that  we  measure 
u  without  errors;  that  is,  the  measurement  of  the  windspeed  u  is  a  true  windspeed.  Obviously 
u  is  not  perfectly  known  without  errors,  so  OLR  gives  a  value  of  windspeed  exponent  7  that  is 
not  statistically  correct.  If  one  were  to  regress  u  vs  o°  instead,  the  value  of  7  obtained  would 
be  larger  than  usually  reported.  Detailed  discussion  of  OLR  can  be  found  in  mathematical 
statistics  books  [20],  [21],  [22]. 

In  contrast.  University  of  Kansas  researchers  often  use  orthogonal  regression  (OR),  rather  than 
OLR.  In  orthogonal  regression,  the  best-fit  line  is  that  for  which  the  mean-square  distance 
normal  to  the  regression  line  is  a  minimum.  The  slope  of  the  best  fit  line,  which  represent  the 
wind  speed  exponent  7,  lies  between  the  low  value  characteristic  of  OLR  for  o°vsu  and  the 
high  value  of  OLR  for  u  vs  o“.  With  weighting  to  equalize  the  standard  deviations,  the  scale 
factor  f  is  simply  the  ratio  of  the  standard  deviations  in  the  two  directions,  and  it  is  also  the 
geometric  mean  of  the  slopes  obtained  by  regression  o°  vs  u,  and  u  vs  o°.  Thus  the  orthogonal 
regression  line  with  slope  b  (7)  must  always  lie  between  the  other  two  regression  lines. 
[23](R.K.Moore,1992). 

The  orthogonal-regression  expression  for  the  slope  b  of  the  regression  line  can  be  shown  to  be: 
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**-2.(0  refers  to  the  point  with  probability  t/2  to  the  right  of  the  Student’s  distribution  function 
with  N-2  degrees  of  freedom.  [22]  (Cramer.  1964). 


For  a  two-dimensional  distribution,  we  write 
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In  particular,  the  quantity  t  defined  by  the  relation 


T 


(12) 


is  the  correlation  coefficient  of  the  sample. 
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In  (25](Boggs,  P.T.  ,1990), [24],  [23]  and  [22],  OR  is  discussed  in  more  detail. 


Another  approach  to  obtaining  the  best  slope  for  a  regression  line  is  due  to  Wald[26](Wald,A. 
1940).  Here  one  divides  the  X  axis  (log  u  in  our  case)  into  two  regions.  One  obtains  mean 
values  of  the  data  for  X  and  Y(o°  dB  in  our  case)  directions  in  each  region.  Then  one  joins  the 
means  with  a  straight  line.  Assuming  the  total  number  of  sample  points  N  is  even,  the  estimator 
in  its  simplest  form  is 


b 


(13) 


[26]  (Wald,  A.,  1940) 

To  compare  the  above  methods,  we  find: 

*  If  there  are  only  measurement  errors  in  one  (Y)  of  the  variables,  the  ordinary  least-square 
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estimation,  OLR,  produces  the  best  results. 

*  The  OR  provides  the  best  estimate  when  x  and  y  have  similar  scales. 

*  For  many  cases  where  measurement  error  occur  in  both  X  and  Y,  the  Wald  method  gives  better 
results  than  either  OLR  or  OR. 

2-2,  Measurement  Errors 

Errors  occur  in  both  radar  and  windspeed  measurements.  Radar  errors  may  be  due  to  noise, 
calibration  problems,  and  errors  in  measuring  the  pointing  angle.  Windspeed  errors  may  result 
from  several  different  causes. 

The  wind  that  affects  the  surface  ripple  structure,  and  therefore  the  radar  signal,  is  at  the  air-sea 
boundary.  However,  wind  is  usually  measured  with  anemometers  placed  some  distance  above  the 
surface.  Standard  heights  for  ships,  for  example,  are  10  m  and  19  m.  When  the  atmosphere  is 
turbulent,  these  measured  winds  are  similar  to  those  at  the  surface.  For  stratified  atmospheres, 
however,  the  surface  winds  can  be  much  less  than  at  anemometer  height.  Some  investigators  use 
model  functions  to  convert  the  anemometer  measurements  to  friction  velocity  u*.  Other 
investigators  convert  measured  winds  at  any  height  to  equivalent  winds  at  10  m  in  a  neutral 
atmosphere.  These  models  depend  on  knowledge  of  the  degree  of  stratification,  usually  obtained 
from  air-sea  temper*:ure  differences.  They  are  far  from  perfect.  See  Figure.2-1.  [27](R.A. 
Vaughan,  1988). 
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Figure  2-1  Linked  hodographs  of  velocity  vectors  in  oceanic  and 

atmospheric  PBLs. 


The  radar  experimental  results  are  given  in  terms  of  either  (1)  measured  winds  at  anemometer 
height,  (2)  equivalent  10-m  neutral-stability-atmosphere  winds,  or  (3)  u*  obtained  by  model 
functions.  For  this  reason,  one  expects  that  the  values  obtained  for  both  A  and  y,  which  are 
magnitude  coefficient  and  wind  speed  exponent,  in 


oB  ■  A  ur 


vary  considerably  between  reports  of  different  experiments. 

Also,  aircraft  cannot  measure  directly  at  the  anemometer  location.  Clearly,  they  must  measure 
o°  over  some  region  that  may  extend  several  km  from  the  anemometer.  The  winds  at  the  aircraft 
measurement  sites  will  be  different,  therefore,  from  those  at  the  anemometer  because  of  small- 
scale  variations  in  the  wind  field. 

For  all  of  those  reasons,  one  expects  that  most  measurements  of  0°  are  less  noisy  than  the  wind 
measurements.  However,  most  investigators  still  use  the  form  of  OLR  that  assumes  perfect  wind 
measurements,  to  calculate  the  wind  exponent  y.  Of  course,  more  meaningful  results  would 
occur  with  OLR  if  the  o°  were  chosen  as  the  independent  value.  This  would  be  true  both 
because  the  relative  error  in  o°  may  be  less  and  because  one  wishes  to  use  o°  to  estimate  u,  not 
u  to  estimate  o°. 
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CHAPTER  THREE 


WINDSPEED  EXPONENT  MEASUREMENT  IN  SAXON-FPN 


SAXON-FPN  was  conducted  at  the  German  Forschungsplatform  Nordsee  in  November,  1990. 
The  high  sea  state  of  the  North  Sea  during  November  provides  excellent  conditions  for 
investigations.  Researchers  of  the  University  of  Kansas  used  the  35-GHz  (Ka-band)  radar  Vector 
Slope  Gauge  (VSG)  and  a  dual-frequency  FM-CW  scatterometer  radar,  operating  at  C  and  X 
bands,  to  make  instantaneous  measurements  of  the  two  orthogonal  components  of  the  slope  vector 
(  along  and  perpendicular  to  the  radar  look  direction).  Unfortunately  ,  the  radar  tracker  of  the 
dual-frequency  scatterometer  did  not  function  properly  under  most  conditions.  Hence,  most  of 
the  data  collected  with  the  X-  and  C-band  radar  could  not  be  used  for  further  analysis.  Also, 
since  an  amplifier  failed  in  the  VSG,  the  data  collected  after  November  27  are  less  useful  for  our 
research  because  of  low  signal-to-noise  ratio.  For  these  reasons,  only  Ka-band  data  (VSG) 
before  November  27  are  discussed  here. 

3-1.  Backscattering  Power  Measurement 

During  the  SAXON-FPN  experiment,  KU’s  systems  were  installed  on  the  northwestern  comer 
of  the  platform  at  about  22  m  altitude.  See  Fig.  3-1.  The  VSG  used  switchable  antenna  beams 
that  illuminated  three  closely  spaced  footprints.  The  three  beams  were  formed  by  a  single 
parabolic  reflector  antenna  with  switched  feeds.  The  antenna's  3-dB  beam  width  was  roughly 
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2®.  and  the  beam-switching  frequency  was  30  Hz.  Fig.  3-2  shows  the  footprint  of  the  VSC.  At 
each  switching  instant,  one  measurement  of  the  range,  and  consequently  the  ocean  wave-height, 
was  obtained.  The  instantaneous  wave  height  at  the  center  of  the  three-beam  footprint  was 
obtained  by  averaging  the  height  measurements  from  all  three  beams.  The  instantaneous 
backscattered  power  was  detected  by  a  mean-square  detector  for  each  beam.  For  more 
information,  see  [28](Bita  B.Sistani,  1993). 

Since  every  run  of  the  radar  observation  was  50  minutes,  over  30,000  data  were  obtained.  Thus, 
in  each  10-minute  interval,  the  mean  power  estimate  used  about  6,000  data.  Because  the  beams 
were  not  identical,  normalization  to  beam  2  requires  multiplying  the  first  and  the  3rd  beams  by 
1.79  and  1.19.  The  normalized  powers  in  each  of  the  three  beams  during  10  minutes  were 
averaged.  After  that,  the  three  averaged- power  values  were  averaged  again  to  obtain  the  mean 
power  of  all  the  three  beams  in  the  10  minutes.  Some  of  these  data  are  shown  in  the  Table  1  and 
Table  2. 
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Table  I 

Part  of  SAXON  Upwind  data 


RUN  No. 

MEAN- 

POWER 

WIND- 

SPEED 

WIND- 

DIR 

RLD-DEG 

WAVE- 

DIR 

INC-ANG 

DEG 

11191047 

1.507 

9.04 

315 

315 

282 

47 

1.409 

8.65 

318 

315 

282 

47 

1.404 

8.32 

314 

315 

282 

47 

1.357 

8.98 

309 

315 

282 

47 

1.506 

8.25 

311 

315 

282 

47 

11191205 

1.231 

7.27 

302 

315 

282 

47 

1.252 

7.65 

294 

315 

282 

47 

1.077 

9.41 

289 

315 

282 

47 

1.318 

8.73 

234 

315 

282 

47 

1.536 

8.73 

283 

315 

286 

47 

11191749 

1.208 

9.24 

272 

315 

284 

47 

1.558 

7.78 

273 

315 

286 

47 

1.583 

8.91 

291 

315 

286 

47 

1.339 

8.05 

301 

315 

286 

47 

1.465 

7.18 

297 

315 

286 

47 
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J  11191846 

1.629 

9.44 

301 

315 

288 

47 

0.977 

6.27 

302 

315 

288 

47 

0.721 

4.67 

287 

315 

288 

47 

0.584 

4.11 

283 

315 

288 

47 

0.340 

4.11 

273 

315 

288 

47 

11192246 

0.885 

6.48 

297 

315 

278 

47 

0.638 

4.99 

300 

315 

278 

47 

0.667 

4.27 

295 

315 

278 

47 

0.467 

4.27 

287 

315 

278 

47 

0.432 

4.11 

296 

315 

278 

47 

11211114 

0.311 

5.80 

58 

45 

356 

45 

0.208 

5.33 

61 

45 

356 

45 

0.222 

5.56 

64 

45 

356 

45 

0.159 

5.56 

69 

45 

351 

45 

0.143 

5.56 

75 

45 

351 

45 

11211213 

0.086 

3.38 

83 

45 

354 

45 

0.094 

2.58 

71 

45 

354 

45 

0.103 

3.07 

59 

45 

354 

45 

0.107 

3.54 

58 

45 

354 

45 

0.190 

4.11 

53 

45 

0 

45 

Table  II 


PART  OF  SAXON  CROSSWIND  DATA 


S  RUN  NO. 

MEAN- 

POWER 

WIND- 

SPEED 

W1ND- 

DIR 

RLD-DEG 

WAVE- 

DIR 

INC-ANG  J 
DEG  | 

11200925 

0.117 

3.95 

55 

315 

279 

47 

0.086 

4.19 

50 

315 

279 

47 

0.084 

4.19 

48 

315 

279 

47 

0.094 

4.19 

53 

315 

279 

47 

0.086 

3.95 

64 

315 

279 

47 

11201733 

0.14 

7.73 

48 

315 

305 

45 

0.151 

7.94 

48 

315 

307 

45 

0.157 

7.37 

46 

315 

307 

45 

0.189 

7.37 

44 

315 

307 

45 

0.169 

7.37 

46 

315 

307 

45 

11201839 

0.194 

7.43 

54 

315 

307 

45 

0.188 

7.10 

52 

315 

315 

45 

0.193 

6.85 

49 

315 

315 

45 

0.163 

7.11 

46 

315 

315 

45 

0.144 

6.62 

46 

315 

315 

45 
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11202305 

0.134 

5.68 

46 

315 

59 

45 

0.124 

5.76 

47 

315 

59 

45 

0.114 

5.68 

41 

315 

53 

45 

0.180 

6.76 

39 

315 

53 

45 

0.133 

6.27 

38 

315 

53 

45 

11230953 

0.480 

8.52 

110 

22 

128 

45 

0.533 

8.71 

107 

22 

128 

45 

0.608 

8.65 

110 

22 

128 

45 

0.586 

8.19 

115 

22 

128 

45 

0.536 

8.98 

120 

22 

128 

45 

11231107 

0.450 

7.35 

124 

22 

128 

45 

0.471 

7.73 

113 

22 

128 

45 

0.396 

7.45 

111 

22 

128 

45 

0.417 

6.86 

113 

22 

128 

45 

0.408 

7.00 

113 

22 

128 

45 

11231243 

0.232 

7.00 

96 

22 

119 

45 

0.245 

6.78 

93 

22 

119 

45 

0.194 

6.72 

93 

22 

119 

45 

0.195 

7.00 

88 

22 

119 

45 

0.203 

7.00 

90 

22 

119 

45 
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3-2.  Wind  Data  Classification 


In  microwave  remote  sensing,  the  ocean  surface  wind  direction  and  speed  as  well  as  the  wave 
direction  are  very  important.  Most  of  the  author's  research  work  in  the  past  two  years  was  to 
study  the  behavior  of  the  radar  scattering  cross-section  under  different  wave  and  wind  conditions. 


Fig  3-3  illustrates  regions  of  different  wind  or  wave  classification  with  an  explanation  of  what 
each  abbreviation  means.  For  example,  when  the  radar  look  direction  (RLD)  has  -t-20°  difference 
from  the  wind  or  wave  direction,  the  radar  is  defined  to  be  looking  Up  wind  or  Up  wave.  To 
help  understand  the  effect  of  wind  and  wave  on  radar  cross-section  0s,  two  cases  for  each  of  the 
following  wind  and  wave  categories  are  discussed  next: 

(i)  Upwind  under  both  Up  and  Cross-wave  or  only  Upwave 

(ii)  Crosswind  under  both  Up  and  Cross  wave  or  only  Upwave 

Since  downwind  and  downwave  data  are  corrupted  by  the  shadow  of  the  tower,  results  from  CD 
(Crosswind-Downwave),  DC  and  Dn  are  not  reported  here.  Also,  for  each  case  mentioned  above, 
only  sample  results  are  presented.  Figure  3-4  shows  the  range  of  wind  conditions  encountered 
from  19  to  23  November. 

(i)  Case  I:  Upwind  tinder  both  Up  and  Cross  wave  or  only  Upwave 
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In  this  case  two  situations  should  be  considered,  (a)  First  we  consider  the  effect  cm  radar  cross- 
section  o°  under  both  Up  and  Cr  wave  conditions.  In  this  situation,  the  values  of  the  wind 
exponent  y  (slope  of  the  linear  equations)  obtained  by  OLR,  OR  and  log  u  vs  o®  regression  are 
higher  than  most  previous  investigators  found  at  lower  frequencies.  Figure  3-5  shows  the  y  was 
3.11  by  OR,  or  2.33  by  OLR  when  wind  speed  varied  from  3.5  to  10  m/s.  The  standard 
deviation  is  only  0.851. 

(b)  Now  consider  the  results  if  the  data  with  Crosswave  are  removed  from  above  situation.  Figure 
3-6  shows  that  the  value  of  y  obtained  by  OLR  or  OR  are  much  lower  than  the  data  that  include 
Crosswave.  The  slopes  of  the  lines  are  1.64  by  OR  and  1.54  by  OLR  while  the  standard 
deviation  is  0.955.  These  slopes  are  close  the  value  1.6  to  2.0  which  most  of  other  investigators 
obtained.{4]  (Moore,  R.K.  vol.  iii.1986). 

(ii)  Case  II:  Crosswind  under  both  Up  and  Crosswave  or  only  Upwave 
This  category  refers  to  the  case  when  the  surface  wind  is  considered  to  be  orthogonal  to  the 
RLD,  while  the  long  wave  is  either  Up  or  Cross.  That  is,  wind  and  wave  directions  are  nearly 
orthogonal  or  close  together,  (a)  If  we  do  not  remove  the  data  with  Cross  wave  from  our 
calculation,  we  obtain  the  higher  values  of  7  shown  in  Figure  3-7.  The  slopes  of  the  lines  are 
2.04  by  OLR  and  3.26  by  OR.  The  standard  deviation  is  0.769.  (b)  However,  when  only  Up 
wave  data  is  considered,  the  resulting  lower  7  and  higher  standard  diviation  are  shown  in 
Figure  3-8.  The  slopes  are  1.004  by  OLR  and  1.161  by  OLR  and  the  standard  deviation  is  0.881. 
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Figure  3-5  Case  I:  (a)  Upwind  under  both  Up 
Cross-wave. 


3-6  Case  I:  (b)  Upwind  under  only  Upwave. 


3-3.  Results 


From  above  discussion,  we  find  the  following: 

i.  Most  measurements  of  0"  are  less  noisy  than  the  wind  measurements.  In  spite  of  this,  most 
investigators  use  OLR.  which  assumes  perfect  wind  measurements.  Thus,  the  results  from  OLR 
can  be  significantly  in  error  for  scattered  measurement  points. 

ii.  The  wave  direction  caused  different  wind  exponents  y  calculated  by  OR  or  OLR.  Thus,  one 
should  consider  both  wind  and  wave  directions  in  estimating  radar  backscatter  or  inverting  to 
obtain  the  wind  vector. 

iii.  The  wind  exponents  y  in  Crosswind  cases  are  higher  than  in  the  Upwind  case. 

iv.  The  highest  y  occurs  for  the  Crosswind,  Crosswave  condition. 

v.  The  wind  exponents  y  independent  wind  speed  but  dependent  wind  direction. 

Cross-wind  exponent  >  Downwind  >  Upwind. 


27 


( dB)  vs  u  (  m/s) 


(  gp  )  oD  uopoas-ssojo 


Figure  3-7  Case  II:  (a)  Crosswind  under  both  Up 
Cross-wave. 


Figure  3-8  Case  II:  (b)  Crosswind  under  only  Upwave. 


CHAPTER  FOUR 


VARIATION  OF  WIND  EXPONENT  WITH  ANGLE  OF  INCIDENCE 

Experimental  results  are  available  at  several  frequencies  and  under  a  variety  of  conditions.  In 
every  case,  the  experimenters  report  their  results  in  terms  of  growth  of  scattering  coefficients 
proportional  to  windspeed  raised  to  some  exponent  y,  with  the  exponent  being  different  for 
different  radar  look  directions  relative  to  upwind.  The  results  from  different  experiments  give 
different  values  of  y.  but  combining  them  permits  us  to  perform  regression  to  get  the  best 
estimate  possible  from  the  data.  Values  of  y  from  several  previous  experiments  are  plotted  in 
Figures  4.1,  in  which  the  reported  values  of  the  windspeed  exponent  vary  quite  widely.  Note  that 
the  values  from  the  KU  Noordwijk  experiment  are  higher  than  those  from  other  experiments 
because  OR  gives  higher  values  than  the  OLR  used  by  other  investigators. 

4-1.  Derivation  of  Bilinear  Regression 

When  one  examines  the  data  in  the  figures,  one  can  reach  the  conclusion  that  the  variation  with 
angle  of  incidence  falls  into  two  regimes:  below  40°  and  above  40°.  Below  40°  the  exponent 
increases  rapidly  with  incidence  angle.  Above  40°  it  increases  more  slowly.  For  this  ieason,  we 
elected  to  use  a  bilinear  regression  technique  to  fit  the  angular  variation  of  y.  Although  there 
is  no  particular  reason  to  expect  the  variation  to  be  linear  in  either  angular  regime,  the  spread  in 
the  data  is  such  that  use  of  a  nonlinear  form  is  not  justifiable. 


Let  us  separate  all  measured  ys  into  two  groups:  the  incidence  angles  are  S  40°  in  a  group  (A) 
and  >  40°  in  group  (B).  Suppose  the  total  number  of  data  is  N.  the  number  of  groups  for  (A) 
is  N,  and  the  number  of  groups  for  (B)  is  N2.  Also,  let  "x“  represent  the  incidence  angle  and  "y" 
the  value  of  the  wind  exponent  y,  we  join  the  two  linear  regression  equations  at  the  point  (X^, 
Y^,  where  ^=40°.  Therefore  the  two  equations  can  be  written  as: 


Y,‘Yt  +  (X,  -Xk)bt  ,  (X,  s  X*) 


(15) 


and 

Y,  «  Yk  *  (X,  -  X*)  b2  ,  ( X ,  fc  X*) 


where,  b,,  b2  are  the  slopes  of  the  regression  lines. 


Now  set. 


s  -  Ed',  -  Y? 


i.l 


(17) 
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If  we  want  to  get  the  minimum  of  S,  we  take  the  partial  derivative  with  respect  to  S. 
Letting, 


as  /  3Yb  =  0 


the  equation  yields 


H  *1  ^ 

Y,Yi~Y>N  ~b i  E*<  -  bi  E xt  *  bx  xk *  bi  V*a  -  o 


Also  setting  the  partial  derivative  of  S  with  b,  equal  to  zero. 


as  /  ab,  =  o 


Then, 


(18) 
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Y,N*b, 


(21) 


E*.  - 

\  > 


+  *i  E*?  2x*Ex. +  Nixl 


E  -  **E  Y> 


(22) 


and, 


Solving  the  above  equation  set,  Y„,  b!  and  b,  can  be  obtained. 


Using  the  data  of  Seasat,  KU  Noordwijk,  Skylab,  etc.  experiments,  the  relationships  with  the 
incidence  angles  for  upwind,  and  downwind  in  Ku  band  for  W  polarization  show  in  Figure  4-1 
and  Figure  4-2,  while  Figure  4-3  and  Figure  4-4  show  the  relationships  for  HH  polarization.  All 
the  slopes  and  the  standard  deviations  for  each  curve  region  are  summarized  in  Table  3 
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Figure  4-1  Bilinear  regression  for  Upwind  under  VV 
polarization. 
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Bilinear  regression  for  Upwind  under  HH 
polarization. 


Table  III 


Bilinear  parameter 


INCIDENT  ANGLE 

0°-- 40° 

408— 70° 

~ 1 

UPWIND-W 

y,=0.040x+0.33 

y2=0.006x+1.71 

mam 

DOWNWIND-VV 

y,=0.055x-0.13 

y2=0.011x+1.64 

Y„  =  2.09 

UPWIND-HH 

y,=0.043x+0.07 

y2=0.015x+1.39 

mam 

DOWNWIND-HH 

y,=0.070x-0.51 

y2=0.021x+1.45 

Yk  =  2.28  j 

4-2.  Discussion 

The  experiments  summarized  here  were  performed  under  different  conditions.  Hence,  it  is 
necessary  to  fine  optimum  values  of  y  at  different  incident  angles  and,  through  the  wind 
exponent  y.  to  obtain  the  ocean  surface  wind.  The  bilinear  regression  plots  in  Fig.4.1  show,  the 
following: 

(i)  The  wind  exponent  y  increases  as  incidence  angle  increases  for  both  upwind  and  downwind. 

(ii)  Wind  exponent  values  at  Ku  band  for  HH  polarization  are  larger  than  for  W  polarization 
for  incidence  angles  0  >  20°.  The  difference  between  them  increases  with  incidence  angle. 
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(iii)  The  difference  of  upwind  and  downwind  is  obvious  after  40s. 


(iv)  The  values  from  the  Seasat  experiment  (14.65  GHz)  and  the  KU  Noordwijk  experiment 
(15.0  GHz),  which  even  are  higher  than  others  because  KU  researchers  used  the  OR  method, 
have  very  similar  values  in  the  bilinear  regression  equation. 
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CHAPTER  FIVE 


CONCLUSIONS  AND  SUGGESTION 


Wind  exponents  (y)  are  very  important  in  the  study  of  radar  scattering  cross-section  under  ocean- 
surface  conditions.  Most  experimenters  apply  the  ordinary  least-square  regression  (OLR). 
However,  OLR  accounts  for  error  in  only  one  variable  in  the  regression,  and  therefore  the  results 
of  OLR  may  be  not  statistically  correct  In  contrast,  the  orthogonal  regression  (OR)  provides  a 
good  way  to  estimate  two  correlation  variables  (x.y)  when  both  of  them  have  errors.  The  best 
results  are  obtained  if  x  and  y  have  similar  scales. 

Our  study  shows  that  OR  has  agreement  with  OLR  when  x  is  perfectly  known  without  errors. 
Otherwise,  OR  and  OLR  have  different  results.  For  example,  SAXON-FPN  experiment  (35 
GHz)  data  demonstrates  that  orthogonal  regression  (OR)  and  ordinary  linear  regression  (OLR) 
methods  have  similar  wind  exponent,  y,  if  the  wave  and  the  wind  directions  are  same.  If 
however,  the  directions  are  not  the  same,  the  estimates  disagree.  In  this  case,  the  estimate  of  y 
from  OR  is  higher  than  OLR.  It  is  also  seen  that  y  for  cross  wind  is  higher  than  y  for  up  wind 
in  both  OR  and  OLR  methods.  Thus,  it  is  seen  that  the  wind  exponent,  y,  has  a  highest  value 
when  the  wind  direction  and  the  wave  direction  are  normal  to  the  radar. 

As  an  important  result,  when  wave  direction  is  similar  with  wind  direction,  the  slop  of  the  best 
fit  line  either  in  OR  or  OLR  is  quite  different  than  that  while  wave  direction  is  orthogonal  wind 
direction.  Therefore,  it  is  obviously  to  see  that  the  wavedirection  is  an  important  factor  which 
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influence  the  wind  exponent  y.  These  above  results  may  explain  some  of  the  scatter  of  y  values 
reported  in  previous  experiments.  Unfortunately,  in  most  cases,  no  wave-direction  measurements 
were  even  made. 

Since  the  1950s,  many  experiments  on  scattering  from  the  ocean  have  been  performed.  They 
give  similar  but  different  results  for  the  wind  exponent  y  that  depends  on  angle  of  incidence. 
Wind  exponent,  y,  appear  to  vary  with  angle  of  incidence  in  two  regions.  The  data  is  varied 
enough  that  they  may  be  represented  by  two  linear  fits.  The  method  of  bilinear  regression  is 
a  simple  but  obvious  way  to  show  the  relation  between  y  and  angle  of  incidence.  From  the  plots 
shown  in  Fig.4-1  and  Fig.4-2,  we  can  have  an  idea  as  to  which  model  function  must  be  used  in 
future  experiments,  such  as  EOS.  The  disadvantage  of  bilinear  regression  is  that  we  have  to 
predetermine  the  vertex  point,  which  can  be  estimated  based  only  on  experience,  and  there  is  no 
particular  reason  to  expect  the  relation  between  y  and  angle  of  incidence  to  be  linear. 

As  bilinear  regression  important  results,  the  y  values  by  OR  are  very  identical  the  bilinear 
regression  equation,  such  as  Seasat  data,  and  the  different  slops  of  upwind  and  downwind  in  the 
bilinear  equations  are  obviously  after  40°.  Also,  the  wind  exponent  y  increases  as  incidence 
angle  increases  for  both  upwind  and  downwind,  but  the  slopes  for  each  region  are  quite  different 

As  a  suggestion  for  future  studies,  if  enough  data  are  obtained,  we  may  use  the  OR  method 
instead  of  the  OLR  method  to  specify  y  exactly. 
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